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ABSTRACT 

A  boundary  collocation  method  is  presented  for  the  computation  of 
eigenvalues  and  eigenfunctions  of  the  Laplacian  on  bounded  simply  connected 
plane  regions  with  smooth  boundary.  Particular  emphasis  is  placed  on  the 
selection  of  an  appropriate  function  space  and  its  basis  for  the  approximation 

of  the  eigenfunctions.  Numerical  examples  are  presented. 

v. 


AMS  (MOS)  Subject  Classifications s  65N25,  65N35 

Key  Words:  Eigenvalue ,  eigenfunction,  Laplacian,  boundary  collocation 
Work  Unit  Number  3  (Numerical  Analysis) 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


SIGNIFICANCE  AND  EXPLANATION 


Several  successful  applications  of  the  boundary  collocation  method  for 
computing  eigenvalues  and  eigenfunctions  of  the  Laplacian  on  certain  regions 
exist  in  the  literature  (see  e.g.  Fox,  Henrici  and  Moler  [5],  and  Moler 
[6]).  However,  Bates  and  Ng  [t]  reported  that  the  method  did  not  converge  for 
some  regions.  We  suggest  a  strategy  for  selecting  collocation  points, 
function  space  for  the  approximation  of  eigenfunctions,  and  basis  of  this 
function  space.  Our  results  apply  to  regions  for  which  there  are  convergence 
difficulties  with  the  methods  referred  to  above. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


A  COLLOCATION  METHOD  FOR  COMPUTING  EIGENVALUES 
AND  EIGENFUNCTIONS  OF  THE  LAPLACIAN 

Lothar  Reichel 

1.  INTRODUCTION 

Let  ft  be  a  bounded  simply  connected  plane  open  region  with  smooth 
boundary  3ft ,  and  let  ft  :=  ft  U  9ft.  We  consider  the  computation  of 
eigenvalues  X  and  eigenfunctions  u  of 

(Au  +  Xu  =  0  in  ft 

u  -  0  on  9ft  * 

Several  successful  applications  of  the  boundary  collocation  method  to  solve 
this  problem  on  certain  regions  exist  in  the  literature,  see  e.g*  Fox,  Henrici 
and  Moler  [5],  and  Moler  [6].  The  purpose  of  this  paper  is  to  describe  a 
strategy  for  choosing  collocation  points  and  function  space  for  the 
approximation  of  eigenfunctions.  We  also  discuss  the  selection  of  basis  of 
the  function  space*  The  need  for  such  a  strategy  is  illustrated  by  the 
computational  results  of  Bates  and  Ng  [ 1 ] .  Choosing  function  space  and 
collocation  points  quite  arbitrarily,  they  not  seldom  experience  convergence 
difficulties  and  conclude  [1],  p.  153,  'The  point  matching  or  collocation 
method  is  an  attractive  economique  technique  from  the  point  of  view  of 
programming  and  computer  time  [••••].  Unfortunately  the  method  is  not  always 
valid*,  ftonconvex  regions  appear  to  be  especially  difficult. 

We  will  represent  eigenf unctions  using  an  integral  operator  introduced  by 
Vekua  [9].  This  operator  maps  functions  analytic  in  ft  onto  solutions  of 
(1.1).  The  connection  between  eigenfunctions  and  analytic  functions,  to  be 
described  in  detail  in  Section  2,  implies  that  the  approximation  of  an 
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eigenfunction  is  equivalent  to  the  approximation  of  a  certain  analytic 
function  on  ft.  The  selection  of  function  space  for  the  approximation  of 
eigenfunctions  can  be  studied  by  considering  the  approximation  of  analytic 
functions  in  ft.  Approximation  of  analytic  functions  in  ft  by  rational 
functions  with  fixed  poles  is  discussed  in  Section  3*  When  approximating 
analytic  functions  by  rational  functions  with  fixed  poles  using  interpolation, 
the  poles  and  interpolations  points  (-  collocation  points)  should  not  be 
allocated  independently.  The  connection  between  analytic  functions  and 
solutions  of  (1.1),  as  well  as  the  reported  difficulties  in  Bates  and  Ng  [1], 
suggest  that  the  choice  of  function  space  for  the  approximation  of 
eigenfunctions  and  the  selection  of  collocation  points  should  not  be  done 
independently  either.  We  begin  by  allocating  collocation  points  on  3ft,  and 
then  define  a  function  space  which  depends  on  the  distribution  of  collocation 
points.  Our  collocations  scheme  is  presented  in  Section  4,  which  also 
contains  a  brief  discussion  on  the  computation  of  error  bounds  derived  by 
Moler  and  Payne  [7).  Numerical  examples  are  found  in  Section  5. 


2.  AN  INTEGRAL  REPRESENTATION 


This  section  summarizes  some  relevant  results  on  the  integral  operator  of 
Vekua,  and  on  the  smoothness  of  the  eigenfunctions.  Let  ft  be  as  in  Section 
1 ,  and  let  denote  the  set  of  functions  whose  pth  derivative 

is  Holder  continuous  on  ft  with  Holder  constant  y .  Let  X  denote  a  real 
nonnegative  constant.  The  following  lemma  is  a  special  case  of  Theorem  4.2  of 
Eisenstat  [3]. 

Lemma  2. 1 

Let  u  e  c2 (£1)  n  p  +  Y  >  0,  satisfy 

(2.1)  Au  +  Xu  ~  0  in  fl  . 

Then  u  has  the  representation 


(2.2)  u(z)  »  Re(4>U)  -  j  U-*t>  /  ♦<C)J1(/ Mz-OU-^Odc),  z  e  fl  , 


where  <M  z)  is  a  function  analytic  in  ft  and  $  e  C^'^tft).  J^(r)  denotes 
the  1st  order  Bessel  function  of  the  1st  kind.  $  is  uniquely  determined  by 
the  requirement  »  0.  z ^  is  an  arbitrary  but  fixed  point  in  ft.  z 

denotes  the  complex  conjugate  of  z ,  and  similarly  for  z^. 

Conversely,  if  $  is  analytic  in  ft  and  $  e  C^'^(ft),  then  u  defined 
by  (2.2)  satisfies  (2.1)  and  u  e  C2(fl)  n  (^'^(5).  ■ 

The  lemma  shows  that  the  approximation  of  an  eigenfunction  of  (1.1)  is 
equivalent  to  the  approximation  of  an  analytic  function  $.  If  3ft  is 
analytic,  then  u  solving  (1.1)  can  be  continued  analytically  across  3ft, 
hence  <p ,  defined  by  lemma  2.1,  can  be  continued  analytically  across  3ft. 

For  nonanalytic  boundary  curves  the  following  lemma  gives  the  smoothness  of 
u  on  3ft.  Hie  lemma  is  a  special  case  of  theorem  7.1  of  Eisenstat  [3]. 

Lemma  2. 2 


Let  3ft  e  C 


max(p, 1 ) ,y 


and  let  u  solve  (1.1).  Then  u  e  CP#^(ft) 


3.  A  METHOD  FOR  APPROXIMATING  ANALYTIC  OR  HARMONIC  FUNCTIONS 

The  connection  between  solutions  of  (1.1)  and  analytic  functions  shown  by 
lemma  2.1  motivates  the  present  section.  We  summarize  an  approximation  method 
for  analytic  functions  that  is  described  in  more  detail  in  [8]  • 

Lemma  3. 1 

Let  T  be  a  smooth  curve  exterior  to  3ft,  and  let  S  denote  the  doubly 
connected  region  between  T  and  3ft.  Let  v  be  a  solution  of 


(3.1) 


Av  =  0  in  S 

<  v  =  1  on  3ft 

v  =  0  on  T 


Then 

(3.2) 


c  :=  /  (  z)  | dz  | 

3ft  3n 


/  ( z) |dz|  >0  , 

r  3n 


where  denotes  the  normal  derivative,  outward  w.r.t.  S.  Let  the  point 

set  {z,  }?  .  be  equidistributed  w.r.t.  c  *  on  3ft,  i.e.  z.  „  is  an 

k,n  k^l  3n  Ln 

arbitrary  point  on  3ft,  and  the  zj  n'  j  *  2,  are  determined  by 


(3.3) 


j  *  1(1)n  , 


where  integration  is  done  counter-clockwise  along  3ft.  Let  the  point  set 

{w,  }"  ]  be  equidistributed  w.r.t.  -c  1  *|^  on  !*•  Let  ^(z)  be  a 
K ,  n  K®  i  d  n 

function  analytic  on  and  interior  to  the  level  curve  :=  {z  s  v(z)  =  w) 

for  some  y ,  0  <  ji  <  1.  Finally  let  r  (z)  be  the  function  in 


(3.4) 


n{ 1,(z  - 


#1,n 


)'\(z 


w  ) 
z  ,n 


'(z  '  wn-1,n)_1} 


,  •  •  • 
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Then 


uniquely  determined  by  interpolating  at  the  nodes  {2^  n^k~1* 

lim  sup  -  r  (z)|  <  e 

n+*  ze3ft  n 

Proof  >  A  proof  based  on  results  of  Walsh  [10]  is  presented  in  [8] •  ■ 

Assume  that  the  nodes  { z.  }?  4  are  given  on  3ft  and  ennumerated  so 

K 1  n  K*“  1 

that  n  precedes  n/  k*1(1)n-1#  when  3ft  is  traversed  counter¬ 
clockwise.  We  wish  to  determine  a  set  of  poles  {w,  }!*  ]  equidistributed 

k  t  n  K*  1 

-1  3v  —  1  3  v 

w.r.t.  -c  *r—  on  a  level  curve  of  v.  Assume  that  c  r  is  a  piecewise 
3n  3n 

constant  function  with  jumps  at  the  The  assumption  that  {z,^  ke 

—1  3v  —1  f  3v  3v 

equidistributed  w.r.t.  c  r—  and  c  J  r~  (z)|dz|  *  1  determines 

3n  aQ  3n  3n 

uniquely  on  3ft.  This  suggests  that  we  determine  level  curves  T  and  point 

sets  {w,  }  ]  by  solving  (3.1)  as  an  initial  value  problem  for  v,  with 
k,n  k=1 

v  and  ~  known  on  3ft.  A  simple  numerical  method  for  solving  this  initial 
on 

value  problem  is  described  in  [8]  and  computed  examples  are  given  in  Section 
5.  Here  we  only  note  that  we  do  not  need  to  solve  the  initial  value  problem 
to  high  accuracy,  and  we  generally  do  not  put  T  very  far  from  3ft. 

Therefore  we  can  obtain  a  satisfactory  approximate  solution  to  this  ill-posed 
problem.  We  next  turn  to  the  choice  of  basis  for  the  representation  of 
rn<z).  The  basis  implicit  in  the  definition  (3.4)  of  Qn  is  generally  ill- 
conditioned.  In  [8]  it  is  shown  that  the  following  basis  functions  are  fairly 
well-conditioned 


v*> 


(3.5)  4  n  z  -  z  n-1  z  -  w 

t  (z)  n  - n  - ks 

k  .  „  z.  -  z,  ,  4  z  -  w . 

k,n  j,n  j»1  j,n 


k  *  1 ( 1 )n  -  1 


This  basis  we  will  use  when  computing  eigenvalues  and  eigenfunctions 


4.  THE  COLLOCATION  METHOD 

The  eigenvalue  problem  (1.1),  we  solve  as  a  sequence  of  Dirichlet 
problems,  i.e.  for  given  X  we  solve 
(4.1a)  Au  +  Xu  =  0 

(4.1b)  u  *  0 

for  u  and  determine  iteratively  values  of  X  for  which  (4.1)  has  nontrivial 

solutions  u.  The  iteration  method  has  previously  been  used  by  Moler  [5] . 

Equations  (4.1)  are  solved  by  boundary  collocation.  Let  {z  and 

k  ,  n  *c—  l 

(w  }**  ]  be  point  sets  on  3ft  and  T,  respectively,  as  described  in 
K ,  n  ks  1 

c  m 

Section  3.  Let  the  set  of  collocation  points  {z,  }.  „ ,  m  =  2n  -  1,  be 

k  ,m  K3  i 

distributed  equidistantly  on  3ft  w.r.t.  to  the  same  density  function  as 

{z  }n  .  Introduce  the  particular  solutions  of  (4.1a), 
k  ,  n  k —  1 


(4.2)  S 


V  (z,X)  s=  Re( 1  -  0.5X (z-z  )  / 
o  i ,  n 


J,</ 

1  ■  at) 


),n  J 


X(z-C) (z-z1 ) 


'2j 


z  J  (/  X(z-C) (z-z  ) ) 

(z,X )  :=  Re(*.(z)  -  0. 5X  (z-z .  )  /  1.(0  - - ;  d?) 

J  ,,n~  J  ' 


1#n 


/ 


X(z-C) (z-z1 ) 


j  =  1 ( 1 )n  -  1  , 


z  J. (/  X(z-C) (z-z.) ) 

V„.  ,(2,X)  :=  Ira( t  (z)  -  0.5XU-Z,  )  /  *.(«)  - - d?) 

O-1  j  •#**  _  J  J 


*1,n 


/ 


X(z-0(z-z) 


j  =  1  ( 1 )  n  -  1  , 


where  the  functions  £^(z)  are  by  (3.5).  We  require  the  boundary 

condition  to  be  satisfied  at  the  collocation  points.  This  yields  the  system 
of  equation 
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2n-1 


(4,3) 


j-1 


a.\\(z,  ,X)  *  0,  k  =  1 ( 1 )m,  m  =  2n  -  1 

3  3  k,m 


We  seek  X  such  that  (4.3)  has  a  nontrivial  solution 

T 

a  =  (a  ,a_ , . . . ,an  , )  ,  or  equivalently,  we  want  the  matrix 
—  12  2n- l 


(4.4)  A  :=  [a.  a.  .  :=  V.(z_C  ,X),  j  =  1(1)2n-1,  k  =  1(1)m,  m  =  2n-1  , 

kj  kj  j  k,m 


to  be  singular.  We  achieve  this  by  computing  the  singular  value  decomposition 
of  A  =  A( X ) , 

A  *  UEVH  , 

where  U,V  are  orthonormal  (2n  -  1)  x  (2n  -  1)  matrices  and 


Z  =  diag(01,02, 


,a„  A,  with  a.  >  o_  >  ...  >  >  0.  We  determine 

2n—l  »  2  zn-i 


X  such  that  a  m  (X)  vanishes,  i.e.  first  we  tabulate 

2n-1  2n-1 

X  o_  AX)  to  find  a  local  minimum,  and  then  by  fitting  quadratics  we 
2n-  i 

2  2 

minimize  a „  (X).  (a_  „(X)  is  a  smooth  function  of  X  while  o  (X) 

2n-1  2n-1  Zn-i 

* 

is  not).  Let  X  denote  a  value  of  X  corresponding  to  a  local  minimum  of 

* 

a  (X).  X  is  an  approximate  eigenvalue  of  (4.1).  A  method  to  estimate 
2n—  1 

* 

or  bound  the  error  in  X  is  discussed  below.  Following  Moler  [5],  we  note 


that  X  satisfies 


min  min  lA(X)al  =  min  lA(X  )al2 
X  loti  =1  ~  lal  =1 


and  that  a  is  the  last  column  of  V.  Hence, 


2n-1 


(4.5) 


u(z)  *  l  a  ^  V  ^  ( z ,  X  ) 


is  an  approximate  eigenfunction.  If  also  a^^tX  ^  is  near  0#  then  X 


7- 


is  an  approximate  double  eigenvalue.  The  coefficients  of  the  other 

eigenfunction  are  in  the  next  to  last  column  of  V. 

The  integrals  in  (4.2),  we  compute  by  integrating  along  3ft  and  using  a 

c  c 

Lobatto  rule  between  each  pair  of  collocation  points  ((z.  )*  (z,..  ))•  We 

K  ,m  K+  1 ,  m 

note  that  the  functions  lAz)  have  to  be  computed  only  once  at  the  Lobatto 

2 

abscissae.  The  minimization  scheme  for  o_  .(X)  is  designed  so  that  we  in 

2n*  i 

* 

each  iteration  know  values  A. ,X  with  A„  <  A  <  X  .  As  the  iterations 

l  r  l  r 

proceed  A^  -  X^  +  0.  This  enables  the  computation  of  the  Bessel  functions  in 
the  integrands  (4.2)  by  linear  interpolation  when  0  <  X  -  X  is  small. 

r  Xr 

Remark  4. 1 

In  order  to  facilitate  the  detection  of  eigenvalues  the  basis  functions  used 

should  be  well-conditioned  in  the  sense  that  >>  0  f°r  ^  between 

* 

local  minima  of  a*  .(A).  Moreover,  for  A  =  A  ,  a  local  minimum  of 
2n—  1 

o  ,<X),  there  should  be  a  k  such  that 
2n-  i 

VX*’  ”  V.(X*>  •  V2(X*’  ■  •••  '  °2„-1(X">  *  0  ■ 

so  that  the  multiplicity  of  the  eigenvalue  can  be  determined  easily.  In  all 
computed  examples  the  basis  (4.2)  satisfied  these  requirements.  ■ 

Remark  4.2 

The  allocation  of  collocation  points  z  can  be  done  quite  arbitrarily.  In 

k  ,m 

all  examples  we  have  distributed  the  collocation  points  equidistantly  w.r.t. 
arc  length.  ■ 

Remark  4.3 

The  choice  of  function  space  for  the  approximation  of  eigenf unctions  is  based 
on  lemma  3.1,  which  discusses  approximation  of  analytic  functions  by  certain 
rational  function.  By  taking  real  and  imaginary  parts,  we  obtain  the  harmonic 
functions  Vj (z,0 ) , • • • ,V2n_^ (z* 0 ) •  Approximation  of  harmonic  functions  by 
boundary  collocation  has  so  far  been  studied  only  when  the  approximating 


functions  are  harmonic  polynomials,  see  Curtiss  [2]*  These  results  suggest 

c  m 

that  the  collocation  points  (z  }  should  be  used  when  approximating  by 

k  ,m  K=  i 

function  (4,2)  with  X  =*  0.  For  X  >  0  no  results  are  known,  and  we  expect 

for  large  values  of  X  least  squares  collocation  to  yield  better 

approximations.  In  least  squares  collocation,  see  Moler  [6],  one  chooses 

m  >  2n  -  1  in  (4.3)  and  solves  (4.3)  in  the  least  squares  sense.  Again  this 

leads  to  the  singular  value  decomposition  of  the  matrix  (4.4),  but  now  one 

2  2  -1 

seeks  local  minima  of  X  ♦  o  1<X)*(//  u  dxdy)  ,  where  u  is  defined  by 

n 

(4.5).  Numerical  experiments  indicate  that,  indeed,  the  approximation  errors 
in  the  eigenvalues  and  eigenf unctions  decrease  when  least  square  collocation 
is  used.  However,  the  decrease  in  error  was  small,  and  in  all  computed 
examples  a  more  efficient  way  to  reduce  the  error  was  to  increase  the  number 
of  basis  functions  and  let  m  =  2n-1.  In  the  numerical  examples  referred  to, 
3ft  was  analytic  and  X  <  10.  ■ 

We  conclude  this  section  with  some  comments  on  the  computation  of  error 

* 

bounds  as  described  in  Moler  and  Payne  [7].  Let  X  be  a  real  nonnegative 

* 

number  and  define  u  (z)  as 


*  n  * 

u  (z)  :=  l  a  V  (z,x  ) 
j-1  3  3 


Let  X  be  the  eigenvalue  of  (1.1)  closest  to  X  ,  and  let  u^(z)  be  the 


eigenfunction  belonging  to  X^.  Moler  and  Payne  [7]  show  how  bounds  for 

*  * 

X  -  X^  and  u  -  u^  can  be  computed.  We  only  give  their  result  on  bounds 

* 

for  X  -  X  * 
k 


Lemma  4. 1 

*  1/2  rtt  *2  \-1 

Let  e  max  |u  (z)|*(area  of  ft)  •  (J  J  (u  )  dxdyj  .  Then  there  is  an 
zedft  ft 

eigenvalue  of  X^  of  (4.1)  in  the  interval 


-9 


/  #  *  2 

The  double  integral  JJ  (u  )  dxdy  we  compute  by  numerical  quadrature, 

ft 

and  we  neglect  round-offs  during  all  computations.  We  therefore  only  compute 
error  estimates*  The  integration  rule  of  the  next  lemma  was  chosen  because  of 
its  simplicity. 

Lemma  4.2 

Let  I  { ( j ,k)  :  {(x,y)  :  |jh  -  x  |  <  |kh  -  yl  <  ^}  C  ft)  .  Then  for  a 
sufficiently  smooth  function  u  which  vanishes  on  dft, 

I //  u2(x,y)dxdy  -  h2  J  u2(jh,kh)|  <  Ch3  , 
ft  ( j ,k)CI 

where  C  is  a  constant. 

Proof .  See  appendix.  ■ 
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5.  NUMERICAL  EXAMPLES 


f  C  ytk 

In  all  examples  the  collocation  points  /.  .  #  m  *  2n-1,  and  the 

)Cf9  JC*  1 

points  { z  „  in  the  definition  of  the  function  (3.5)  are  equidistributed 

k,n  k*1 

on  3ft  w.r.t.  arc  length.  By  solving  (3.1)  as  an  initial  value  problem  for 

3  v  ^ 

v  and  ^  on  3ft  we  compute  approximate  level  curves  V  of  v  and  point 

sets  {w^  n}j^  which  are  approximately  equidistributed  w.r.t.  -c  on 

?.  Below  we  show  graphs  with  ft/  curves  T  and  points  {w  }!*)•  ^  all 

examples  the  integrals  in  (4.2)  have  been  computed  by  integrating  along  3ft 

and  using  a  10-point  Lobatto  rule  between  each  pair  of  collocation  points 

(z°.  ,zC  «  ),  j»1(1)m-1.  The  double  integrals  in  the  error  bound  of 

j  #m  1 

Lemma  4.1  have  been  computed  by  the  method  of  Lemma  4.2.  We  have  computed 

2 

//  u  dxdy  for  different  values  of  h,  and  decreased  h  until  we  were 

ft 

certain  that  the  value  of  the  integral  was  known  with  2  significant  digits. 

All  computations  were  done  on  a  VAX  11/780  in  double  precision/  i.e.  with  12 
significant  digits.  In  the  first  example,  we  determine  eigenvalues  for  an 
elliptic  membrane  for  which  the  lowest  eigenvalue  is  know  with  high  accuracy. 
Ex.  5.1  ft  :*  {(x,y)  :  —  x  +  y  <  1}  .  Figure  5.1  shows  ft/  approximate 
level  curves  V  and  the  points  wj  3/  j  *  1(1)8/  at  the  vertices  of  the 
outermost  approximate  level  curve.  The  wj,3  are  use<*  define  17  basis 
functions  (4.2). 

Figure  5.2  is  similar  to  Figure  5.1,  but  now  16  points 
j  *  1(1)16  are  shown  on  the  outermost  approximate  level  curve  T.  The  16 
points  Wj  17  define  33  basic  functions  (4.2). 

The  €  in  the  table  is  defined  in  lemma  4.1.  The  table  gives  the  3 
smallest  eigenvalues  determined,  the  smallest  of  which  is  known,  see  [5],  to 
lie  in  the  interval  [3.56672658,  3.566726614]. 


No.  of 
Basis  Facts 


Computed 

Eigenvalue 


e 


Est.  Upper  Bound  for  Eigenvalue 
Est.  Lover  Bound  for  Eigenvalue 


17 


33 


3.56672688  8.5  •  10~5 


3.566726603  2.8  •  10~9 


3.5670 

3.5664 

3.56672662 

3.56672659 


17  6.27535  4.7  •  10”3 


33  6.275432466  9.4  •  10~7 


6.306 

6.245 

6.275439 

6.275426 


17 


10.034 


1.2  *  10 


-1 


33 


10.028402  2.2  •  10 


-5 


11.36 

8.98 

10.02863 

10.02817 


The  error  estimates  generally  appear  to  be  quite  conservative. 


Ex.  5.2  Let  ft  be  the  interior  of  the  curve 

3ft  :*  (x(t)  +  iy(t),  x(t)  tm  1.4  cos(t)  +  2.8  cos(2t)  -  2.8  , 
y(t)  :*  1.6  8in(t)  +  2.45  sin(t  -  0.2)  ♦  0.98  sin(2t)  - 
-  0.49  sin(4t )  -  0.74,  0  <  t  <  2t}  , 
see  Figure  5.3.  The  lovest  eigenvalue  has  been  computed  with  different 
numbers  of  basis  functions 


i 


i 
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HO.  Of 
Basis  Pacta 


Computed 

Eigenvalue 


€ 


Eat.  Upper  Bound  for  Eigenvalue 
Eat.  Lower  Bound  for  Eigenvalue 


33  1.2259  2.4 

49  1.22980  1.4 

1.241 

65  1.2305134  8.4  •  10“3  - 

1.220 

1.2353 

97  1.23053147  3.8  •  10~3  - 

1.2252 

The  pointa  wj#n  for  different  baaea  are  8hown  in  Pigurea  5. 3-5. 6.  The 
allocation  of  polea  for  33  and  65  ba8ia  functiona  is  similar,  and  so  ia 

the  allocation  for  49  and  97  basis  functions.  The  example  shows  that  the 
allocation  of  poles  does  influence  the  rate  of  convergence,  but  the  computed 
eigenvalues  are  not  very  sensitive  to  differences  in  allocation.  ■ 
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appendix 


Proof  of  lemma  4.2 

w.  for  .illicit,  ttat  »ec4(5>.  Lot  s  b.  th.  -qu*r. 

1„„)  ,  1,1  <  lyl  <  |1  *«•»  (  -  »2.  Then,  (41.  =h.  « 


Tf 


JJ  f(x,y)dxdy  * 

s 


( 2f(0,0)  +  flh.O)  *  f (-h,0)  +  f(0,h)  +  f(0,-h))  +  0(h6) 


This  yields 


Tf  =.  h2f(0,0l  *  h4(^f  (0.01  *  A  10.01)  *  0(h6) 

%y~  3y 


Let  S  := 


u  (x,y)  *  Ijh  -  x|  <  §,  Ikh  -  gl  <  n>en 


(j,k)ei 


If  f(x,y)dxdy  »  //  f(x,y)dxdy  +  JJ^  f(x,y)dxdy 


fl\S 


h2  l  f<  jh,kh)  +  h4  I  (jh.kh)  +  fj  (jh.kh))  + 

<j,K)ei  3x  *y 


+  //  f(x,y)dxdy  +  0(h  )  • 
ft\S 

2  3  f  3  f  B  jf flSl^  +  uAu)  yields 

He  now  make  use  of  that  f-u.  2  +  a2  ','‘3x'  ^3y 

3x  9y 
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h2  l  (^f  (Jh.Kh)  Oh.*!.))  - 

(j,k)ei  3x  ®Y 

-  //  2i |  (X,y)  +  ^-f  <x,y>  dx«*y  +  °^2)  " 

*  3x2  3y 

s 

- 2  IJ  (f;)2  *  (57)J  <'"Suaxay  * 0<h2>  ‘  2  '•  “  *■  *  *  0lh*’  ’ 

s 

„h„.  i-  danota.  th.  outward  nona.1  darlw.tlw.  fro.  S.  By  definition  of 

i.  t  Jai.t.nc.  tat««  M  .»d  «  U  <  —  «*-  “  l*  b-“  “ 

«  ao  «.  have  if  I2  dsl  <  d,h,  for  some  constant  dr 

ft,  and  u  =  0  on  <*»*,  we  nave  3U  1 

iS 

Si.il.rly  //  u2d,dy  <  d2n2  //_  and,  <  d,h’.  for  conat»t.  d2.d,. 

Q\8  a\s 

This  completes  the  proof. 


LR/ed 


w  j  f. 
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